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(Dated: March 30, 2005) 

We present the first measurement of the Ai and A3 angular coefficients of the W boson produced 
in proton-antiproton collisions. We study W — > eu e and W — > [xv^, candidate events produced in 
association with at least one jet at CDF, during Run la and Run lb of the Tevatron at y/s—1.8 TeV. 
The corresponding integrated luminosity was 110 pb _1 . The jet balances the transverse momentum 
of the W and introduces QCD effects in W boson production. The extraction of the angular 
coefficients is achieved through the direct measurement of the azimuthal angle of the charged lepton 
in the Collins-Soper rest-frame of the W boson. The angular coefficients are measured as a function 
of the transverse momentum of the W boson. The electron, muon, and combined results are in good 
agreement with the Standard Model prediction, up to order in QCD. 

PACS numbers: ll.30.Er, ll.80.Cr, 12.15.-y, 12.38.-t, 12.38.Qk, 13.38.Be, 13.87.-a, 13.88,+c, 14.70. Fm 



I. INTRODUCTION 

Measurements of the W boson differential cross sec- 
tion, as a function of energy and direction, provide in- 
formation about the nature of both the underlying elec- 
troweak interaction, and the effects of chromodynamics 
(QCD). This differential cross section can be expressed 
us a function of the helicity cross sections of the W, al- 
lowing us to study the W polarization and associated 
asymmetries. Because of the difficulties in fully recon- 
structing a W boson in three-dimensions at a hadron 
collider, the complete angular distribution of the W has 
not been determined yet. In this paper we present the 
first measurement of two of the four significant leading 
angular coefficients of the W boson produced at a hadron 
collider. 

The total differential cross section for W boson pro- 
duction in a hadron collider is given by 

da 3 da u 2 

d{p^) 2 dydcos9d ( j ) ~ 16n d(pf) 2 dy [[ + °° S ' 

+ -A (l-3cos 2 9) + Ai sin 29 cos 4> 

+ ^A 2 sin 2 9 cos 2(f) + A 3 sin cos 

+ A A cos 9 + A 5 sin 2 9 sin 2<j) 

+ A e sin 29 sin <f) + A 7 sin 9 sin </>] (1) 

where pY^ and y are the transverse momentum and the 
rapidity of the W in the laboratory frame, and 9 and 
(j) are the polar and azimuthal angles of the charged lep- 
ton from W boson decay in the Collins-Soper (CS) frame 
[1]. The factors A t (p^,y) are the angular coefficients of 
the W boson, which are ratios of the helicity cross sec- 
tions of the W and its total unpolarized cross section 
da u /d{pf) 2 dy. The CS frame [2] is the rest-frame of the 
W with a z-axis that bisects the angle between the pro- 
ton direction and the direction opposite that of the an- 
tiproton (Figure 1), and it is used because in this frame 
we can in principle exactly reconstruct the azimuthal an- 
gle 4> and the polar quantity |cos#|. Our ignorance of 
the W boson longitudinal momentum, which is due to 
our inability to measure the longitudinal momentum of 
the neutrino, only introduces a two-fold ambiguity on the 
sign of cos 9. It is common to integrate Equation (1) over 



y and study the variation of the angular coefficients as a 
function of p^ . 

To study the angular distribution of the W we must 
choose a particular charge for the boson. In this paper 
we consider the W~ bosons; the W + bosons in our sam- 
ples are CP transformed to be treated as W~ bosons. 
The angular coefficients for the W + are obtained by CP 
transforming Equation (1) [3]. 

If the W is produced with no transverse momentum, it 
is polarized along the beam axis, due to the V-A nature of 
the weak interactions and helicity conservation. In that 
case A A is the only non-zero coefficient. If only valence 
quarks contributed to W production, A4 would equal 2, 
and the angular distribution given by Equation (1) would 
be ~ (1 + cos6*) 2 , a result that was first verified by the 
UA1 experiment [4]. 

If the W is produced with non-negligible transverse 
momentum, balanced by the associated production of 
jets, the rest of the angular coefficients are present, and 
the cross section depends on the azimuthal angle <j> as 
well. The last three angular coefficients A 5 , A e , and A-j 
are non-zero only if gluon loops are present in the pro- 
duction of the W boson. Hence, in order to study all the 
angular coefficients and associated helicity cross sections 
of the W boson in a hadron collider, we must consider 
the production of the W with QCD effects up to order 
a 2 . 

The importance of the determination of the W angular 
coefficients is discussed in [5], and summarized here. It 
allows us to measure for the first time the full differential 
cross section of the W and study its polarization, since 
the angular coefficients are directly related to the helicity 
cross sections. It also helps us verify the QCD effects in 
the production of the W up to order a 2 . For example, 
according to the Standard Model (SM), A 2 is not equal 
to Aq only if the effects of gluon loops are taken into 
account. In addition, A3 is only affected by the gluon- 
quark interaction and its measurement can be used to 
constrain the gluon parton distribution functions. More- 
over, the next-to-leading order angular coefficients A 5 , 
Aq, and Ay are P-odd and T-odd and may play an im- 
portant role in direct CP violation effects in W produc- 
tion and decay [6]. Finally, quantitative understanding 
of the W angular distribution could be used to test new 
theoretical models and to facilitate new discoveries. 
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FIG. 1: Transforming from the laboratory frame to the 
Collins-Soper frame. We first boost to the W rest-frame, then 
rotate the x — z plane so that it coincides with the p-p plane. 
Finally we rotate the frame around the y-axis so that the z- 
axis bisects the angle between p and —p. The positive y axis 
is selected to have the same direction as pes X Pes ■ 



In this paper we present the first measurement of the 
A2 and A3 angular coefficients of the W boson. These 
coefficients fully describe the azimuthal differential cross 
section of the W boson, and they are two of the four 
significant coefficients that describe the total differential 
cross section of the W, given that A\ and the next-to- 
leading order angular coefficients have considerably lower 
values [5, 7]. This measurement is accomplished using 
the azimuthal angle of the charged lepton in the CS W 
rest-frame [8] , and is presented as a function of the trans- 
verse momentum of the W boson. The CS polar an- 
gle analysis is more sensitive to the Aq and A4 angular 
coefficients (see [9, 10] for a measurement of Aq). Be- 
cause Equation (1) arises solely from quantum field the- 
ory, without input from any specific theoretical model of 
W boson production, our experimental results are thus 
model-independent. 



II. THE CDF DETECTOR AND EVENT 
SELECTION 

A. The CDF detector 

The CDF detector is described in detail in [11]. It is 
a general purpose detector of charged leptons, hadrons, 
jets, and photons, produced from proton-antiproton col- 
lisions at the Tevatron accelerator at Fermilab. The W 
and Z bosons are detected through their decay leptons, 
while the transverse momentum of the neutrinos is esti- 
mated from the missing transverse energy of the events 
{$T )• 



The z-axis of the detector coincides with the direction 
of the proton beam and defines the polar angle 6\ a \> in the 
laboratory frame. The y-axis points vertically upward 
and the £-axis is in the horizontal plane, so as to form 
a right-handed coordinate system. The pseudorapidity, 
7?iab = — bi[tan(#i a b/2)], and the azimuthal angle 4>\ab are 
used to specify detector physical areas. 

The tracking system of CDF consists of the silicon ver- 
tex detector (SVX), the vertex time projection chamber 
(VTX) and the central tracking chamber (CTC), all im- 
mersed in a 1.4 T magnetic field produced by a super- 
conducting solenoid of length 4.8 m and radius 1.5 m. 
The SVX, a four layer silicon micro-strip vertex detec- 
tor, is located immediately outside the beampipe. It is 
used to find secondary vertices and provides the impact 
parameter of tracks in the transverse r — </>i a b plane. The 
VTX, located outside the SVX, is a vertex time projec- 
tion chamber that provides r — z tracking information 
up to a radius of 22 cm and pseudorapidity |?7i a b| < 3.5. 
It measures the z-position of the primary vertex. Fi- 
nally, surrounding the SVX and the VTX is the CTC, a 
3.2 m long cylindrical drift chamber containing 84 layers 
of sense wires arranged in five superlayers of axial wires 
and four superlayers of stereo wires. The axial super- 
layers have 12 radially separated layers of sense wires, 
parallel to the z-axis, that measure the r — 4>\ab posi- 
tion of the tracks. The stereo superlayers have six layers 
of sense wires with alternate ~ ±3° stereo angles with 
respect to the bcamlinc, and measure a combination of 
t — <^iab an d z information. The stereo and axial data are 
combined to reconstruct the 3-dimcnsional track. The 
CTC covers the pseudorapidity interval |r/i a b| < 1-0 and 
transverse momentum px > 0.4 GeV [12]. The com- 
bined momentum resolution of the tracking system is 
5p T /pT = ^(O.OOOdpr) 2 + (0.0066) 2 , where p T is the 
transverse momentum in GeV. 

The solenoid is surrounded by sampling calorimeters 
used to measure the electromagnetic and hadronic en- 
ergy of electrons, photons, and jets. The calorimeters 
cover the pseudorapidity range |?7iab| < 4.2 and the az- 
imuthal angle range < 0i a b < 2ir. They are segmented 
in Viib — 0iab towers pointing to the nominal interaction 
point at the center of the detector. The tower granu- 
larity is (A?7i a b x A^i a b) = (0.1 x 15°) in the central 
region (0 < |ryi a b| < 1.1) and (0.1 x 5°) in the plug 
(1.1 < \vub\ < 2.4) and forward (2.4 < \r) Uh \ < 4.2) 
regions. Each region has an electromagnetic calorime- 
ter (CEM in the central region, PEM in the plug re- 
gion, and FEM in the forward region) followed by a 
hadron calorimeter at larger radius from the beam (CHA, 
PHA, and FHA respectively). The central calorimeters 
are segmented in 24 wedges per each half of the de- 
tector (-1.1 < 7? lab < and < ?/i ab < 1.1). The 
CEM is an 18 radiation length lead-scintillator stack 
with a position resolution of 2 mm and an energy res- 
olution of SEt/Et = ^{13.5%/VEtV + (2%F> where 
Et is the transverse energy in GeV. Located six radiation 
lengths deep inside the CEM calorimeter (184 cm from 
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the beamline), proportional wire chambers (CES) with 
additional cathode strip read-out provide shower posi- 
tion measurements in the z and r — 0i a b directions. The 
central hadron calorimeter (CHA) is an iron-scintillator 
stack which is 4.5 interaction lengths thick and provides 
energy measurement with a resolution of SEt/Et = 
V(50%/V£t) 2 + (3%) 2 , where E T is the transverse en- 
ergy in GeV. 

The central muon system consists of three components 
and is capable of detecting muons with transverse mo- 
mentum pr > 1.4 GeV and pseudorapidity |r?i a b| < 1-0. 
The Central Muon Chambers (CMU) cover the region 
|f?iab| < 0.6 and consist of four layers of planar drift cham- 
bers outside the hadron calorimeter, allowing the recon- 
struction of the muons which typically pass the five ab- 
sorption lengths of material. Outside the CMU there are 
three additional absorption lengths of material (0.6 m of 
steel) followed by four layers of drift chambers, the Cen- 
tral Muon Upgrade (CMP). The CMP chambers cover 
the same pseudorapidity region as the CMU, and they 
were introduced to limit the background caused from 
punch-through pions. Finally, the Central Muon Exten- 
sion chambers (CMX) cover the region 0.6 < ?7i a b < 1.0. 
These drift chambers are sandwiched between scintilla- 
tors (CSX). Depending on the incident angle, particles 
have to penetrate six to nine absorption lengths of mate- 
rial to be detected in the CMX. The particle candidate 
stub provided by the muon system is matched with a 
track from the CTC in order to successfully reconstruct 
a muon. 



B. The CDF triggers 

CDF has a three-level trigger system designed to select 
events that can contain electrons, muons, jets, and 1£t ■ 
The first two levels are implemented in hardware, while 
the third is a software trigger which uses a version of the 
offline reconstruction software optimized for speed and 
implemented by a CPU farm. 

At level- 1, electrons were selected by the presence of an 
electromagnetic trigger tower with energy above 6 GeV 
(Run la) or 8 GeV (Run lb), where one trigger tower 
consisted of two adjacent physical towers (in pseudora- 
pidity). Muons were selected by the presence of a track 
stub in the CMU or CMX, where there was also signal 
in the CMP. 

At level- 2, electrons satisfied one of several triggers. In 
Run la, the event passed the trigger if the energy clus- 
ter in the CEM was at least 9 GeV with a seed tower 
of at least 7 GeV, and a matching track with px > 9.2 
GeV was found by the Central Fast Tracker (CFT), the 
fast hardware processor that matched CTC tracks in the 
f ~ '/'lab plane with signals in the calorimeters and muon 
chambers. It also passed the trigger if there was an iso- 
lated cluster in the CEM calorimeter of at least 16 GeV. 
The most common Run lb level- 2 electron trigger requires 
the existence of a cluster in the CEM with at least 16 



GeV and the existence of a matching track in the CFT 
with pt > 12 GeV. The muon trigger at level- 2 required 
a track of at least 9 GeV (Run la) or 12 GeV (Run lb) 
that matched a CMX stub (CMX triggers), both CMU 
and CMP stubs (CMUP triggers), or a CMU stub but 
no CMP stub (CMNP triggers). 

At level-3, reconstruction programs performed 3- 
dimensional track reconstruction. In the Run la level- 
3 electron trigger, most of the accepted events passed 
the requirement that the CEM cluster had Et > 18 
GeV, and was associated with a track of pr > 13 
GeV. The transverse energy of the cluster is defined as 
Et = Esind, where E is the total energy deposited in 
the CEM, and 6 is the polar angle measured from the 
event vertex to the centroid of the cluster. Cuts were 
applied on the shape of the electron shower profile and 
the energy deposition patterns. In the Run lb level-3 
electron trigger, CEM E T > 18 GeV and CFT p T > 13 
GeV requirements were applied. The muon trigger at 
level-3 required that the CFT transverse momentum was 
greater than 18 GeV, the energy deposited in the hadron 
calorimeter was less than 6 GeV, the energy deposited 
in the electromagnetic calorimeter was less than 2 GeV, 
and the extrapolated CTC track was no more than 2 cen- 
timeters away from the muon stub in the CMU chambers 
and 5 centimeters in the CMP or CMX chambers in the 
x direction. Events that pass the level-3 trigger were 
recorded to tape for offline analysis. 



C. The datasets 

The events passing the three levels of our trigger sys- 
tem constitute the inclusive high-p^ electron and muon 
data samples. We apply kinematic and lepton identifica- 
tion cuts, described in Sections II C 1 and II C 2 to obtain 
the inclusive W electron and muon datasets respectively. 
Using these datasets we arrive at the W+jet datasets by 
applying the jet selection cuts described in Section II C 3. 



1. Inclusive W Electron Selection 

After passing the three levels of trigger requirements, 
the following event selection cuts are applied to the in- 
clusive electron data sample: 

• The event must belong to a good run. 

• E T >20 GeV, 

where E T is the transverse energy of the CEM cluster, 
corrected for differences in response, non-linearities, and 
time-dependent changes. 

• ICbl < i, 

where 77j e ab is the pseudorapidity of the electron. 

• The electron must fall in a fiducial part of the CEM 
calorimeter. 

• ISO (0.4)= El^ A /E T lustcI < 0.1, 

where ^Afi= s o 4 ^ s t Qe exc ess transverse energy in a cone 
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of size AR — \J (A0i ab ) 2 + (A?7i a b) 2 centered on the di- 
rection of the electromagnetic cluster, and E^ ustcr is the 
transverse energy of that cluster. 

• £ HAD /£ EM < 0.055 + 0.00045£ e , 

where £; HAD is the energy deposited in the hadron 
calorimeter, and E EM is the energy deposited in the elec- 
tromagnetic calorimeter. 

nmeas pcxp 

• ^hr = o.i4 Ei V(0 , 14) ^; (A£rP)2 < 0.2, ^ 

where Lshr is the lateral shower profile, i^ mcas is the 
energy measured in the « th -tower adjacent to the seed 
tower, E^ Kp is the expectation for the energy in that 
tower, Ai^ exp is the uncertainty on the expected energy, 
and 0.14\/ E meaa is the uncertainty in the measurement 
of the cluster energy. 

• Xces < 10- 

We measure the shower profile along the z direction us- 
ing the CES strips and the shower profile along the x 
direction using the CES wires. By comparing the mea- 
sured x-shape and z-shape to the ones determined from 
test-beam studies we extract the chi-squared quantities 
for the two directions. The chi-squared we use is the 
average of the two. 

• 0.5 < E c /p c < 2.0, 

where E c is the corrected energy of the electron, and p c 
is the beam-constrained momentum of the electron, i.e., 
the momentum determined when the fit trajectory of the 
CTC hits is constrained to pass through the beam line. 

• | AX | < 1.5 cm and \AZ\ < 3.0 cm, 

where AX and AZ are the difference in the x and z 
directions respectively, between the extrapolated CTC 
track and the CES position of the shower. 

• I-^vtxI < 60 cm, 

where ZyTx is the z position of the primary vertex. 

• Photon conversions are removed. 
We next apply the following cuts: 

• $ T > 20 GeV, 

where $r is the missing transverse energy in the event, 
calculated from the energy imbalance in the calorimeters, 
with a correction for the unclustered energy - calorime- 
ter energy not taken into account by the jet clustering 
algorithm - and possible presence of muons. 

• M*f > 40 GeV, 

where Mjf is the W transverse mass. This cut removes 
the background from W bosons decaying into tau leptons 
which subsequently decay into electrons. 

• The event must not be consistent with a Z decaying 
into two observed leptons, or a Z in which one of the 
decay tracks has not been identified. 

The 73363 events passing these cuts constitute our in- 
clusive W electron data sample (Run la: 13290 events 
and Run lb: 60073 events), corresponding to an inte- 
grated luminosity of 110 pb _1 (Run la: 19.65±0.71 pb _1 
and Run lb: 90.35 ± 3.70 pb" 1 ). 



2. Inclusive W Muon Selection 

After passing the three levels of trigger requirements, 
the following event selection cuts are applied to the in- 
clusive muon data sample: 

• The event must belong to a good run. 

• p£ > 20 GeV, 

where pj< is the beam-constrained transverse momentum 
of the muon (determined by a fit to the CTC hits, con- 
strained by the beam line). 

• The muon must be fiducial and central (pscudora- 
pidity |< ab | < 1). 

• ISO(0.4)= ^a^oVPt < where e ar=oa is 
th e excess transverse energy in a cone of size AR = 
^/(A0i a b) 2 + (Afyiab) 2 centered on the direction of the 
muon. 

• Shad < 6 GeV, 

where -Ehad is the energy deposited in the hadron 
calorimeter tower traversed by the muon. 

• E EM < 2 GeV, 

where Eem is the energy deposited in the electromagnetic 
calorimeter tower traversed by the muon. 

• |AV C mu| < 2 cm, |AV C mp| < 5 cm, |AV C mx| < 5 
cm, 

where AXcmu , AXqmp and AXcmx are the differences 
between the x position of the stub in the muon chambers 
and the extrapolation of the CTC track to these muon 
chambers. 

• I-^vtxI < 60 cm. 

• The event must pass the cosmic ray filter. 

• The impact parameter must be \do\ < 0.2 cm. 

• \Z - Zvtx\ < 5 cm, 

where Z is the z-position of the muon track. This cut, 
combined with the previous two, significantly reduces the 
cosmic muon background. 

We next apply the following cuts, as in the electron 
case: 

• pr > 20 GeV. 

• > 40 GeV. 

• The event must not be consistent with a Z decaying 
into two observed leptons, or a Z in which one of the 
decay tracks has not been identified. 

The 38601 events passing these cuts constitute our in- 
clusive W muon data sample [Run la (CMUP): 4441 
events, Run la (CMNP): 955 events, Run lb (CMUP): 
20527 events, Run lb (CMNP): 3273 events, and Run lb 
(CMX): 9405 events], corresponding to an integrated lu- 
minosity of 107 pb- 1 [Run la (CMUP): 18.33±0.66 pb" 1 , 
Run la (CMNP): 19.22 ± 0.69 pb" 1 , Run lb (CMUP): 
88.35 ± 3.62 pb" 1 , Run lb (CMNP): 89.20 ± 3.66 pb" 1 
and Run lb (CMX): 88.98 ± 3.65 pb" 1 ]. 



3. Inclusive W+jet Event Selection 

Our final analysis dataset consists of those W events 
which include at least one jet with > 15 GeV, 

1^*1 < 2.4, and AR}^ > 0.7, where AR}^ = 
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TABLE I: The electron (N e ) and muon (iV M ) CDF data event 
yields for inclusive VK+jet production. The muon event yields 
are lower relative to the electron event yields, due to lower 
muon efficiencies and acceptances. 



Data event yields for inclusive VK+jet production 



Pt (GeV) 



15-25 
25-35 
35-65 
65-105 



5166 
3601 
3285 
624 



2821 
1869 
1880 
371 



(Afjji^.) 2 + (A^.) 2 , and AryJ^. and At^. are the 

differences in pseudorapidity and polar angle between the 
charged lepton and the jet in the laboratory frame. The 
results of the analysis pertain to the W~ boson. All W + 
bosons in the sample are CP transformed to be treated 
as W~ bosons [13]. 

These requirements leave 12676 electron W^+jet events 
and 6941 muon VK+jet events, with 15 < p% < 105 
GeV, where pY^ is the transverse momentum of the W 
boson, defined as the vector sum of the IpT and charged 
lepton transverse momentum. The data event yields for 
the four p% bins (15 < p% < 25 GeV, 25 < p% < 35 
GeV, 35 < pf < 65 GeV, and 65 < p^ < 105 GeV) arc 
presented in Table I. 

The actual number of events is not of critical impor- 
tance for us, because we are interested in the shape of the 
distributions and not the absolute event yields. We thus 
analyze the distributions normalized to unity. We will 
come back to the actual event yields after the inclusion 
of the background (Section V) and systematic uncertain- 
ties (Section IX). 



III. THE MONTE CARLO SIMULATION 

A. The DYRAD Monte Carlo event generator 

DYRAD [14] is the next-to-leading order VK+jet event 
generator used to establish the SM prediction. We in- 
clude the "1-loop" processes, since these affect the next- 
to-leading order angular coefficients and more completely 
simulate the events we study. This generator is of order 
a 2 in QCD, generating up to two jets passing the mini- 
mal requirement of > 10 GeV if the Feynman dia- 
gram does not contain any gluon loops, and generates up 
to one jet with the same requirement if a gluon loop is 
present in the Feynman diagram. As a result, DYRAD 
does not appropriately model events with more than two 
jets. These extra jets in the data occupy low and high 
values of the azimuthal angle <f> in our CS frame. We 
are careful not to bias our measurement due to this ef- 
fect (see Section VIII). The jet transverse energy cut 
of 10 GeV is required because the theoretical calcula- 
tions are unreliable for small jet transverse energies due 
to infrared and collinear divergencies. A jet-jet angular 



separation cut of greater than 0.7 in f?i a b-<^iab space is im- 
posed, which is important for the definition of a jet. No 
additional kinematic cuts for the jet, charged lepton, and 
neutrino are required in order to obtain a reliable theo- 
retical prediction of the angular distribution of the W. 
The cross section for inclusive VF+jet production calcu- 
lated up to order a 2 is 722.51 ± 3.89 pb for W~ and W+ 
bosons combined. This simulation uses Q 2 — (M^? lc ) 2 , 
where M^ lc = 80.3 GeV is the pole mass of the W, 
CTEQ4M(A = 0.3 GeV) parton distribution functions 
[15], and 0.7-cone jets in ryiab-^iab space. 

In order to obtain smooth SM kinematic distributions 
up to p^ — 100 GeV, and especially smooth cos 9 vs. 
(j) distributions for different p^ regions, we generated a 
large sample of DYRAD events (-250 M). This Monte 
Carlo event sample size was required since events with 
negative weights, corresponding to the gluon loop ma- 
trix elements, produce significant fluctuations in the kine- 
matic distributions with limited statistics. The DYRAD 
simulation allows us to establish the SM prediction for 
the (f> distribution of the charged lepton and the predic- 
tions for the angular coefficients and helicity cross sec- 
tions of the W up to order a 2 [7]. The expected <f> dis- 
tributions for four p^ bins are shown in Figure 2. For 
zero p^ we expect a flat distribution, whereas the QCD 
effects at higher p^ result in two minima. In order to 
simulate the detector response, we pass the generator 
events through the fast Monte Carlo detector simulator, 
described in the next section. 



B. The fast Monte Carlo detector simulation 

The fast Monte Carlo (FMC) CDF detector simulation 
includes the detailed geometry of the detector, geometri- 
cal and kinematic acceptances of all subdetcctors, detec- 
tor resolution effects parameterized using gaussians ob- 
tained explicitly from data, detailed magnetic field map, 
and multiple Coulomb scattering effects. The integrated 
luminosities, lepton identification and trigger efficiencies, 
and all experimental cuts imposed on the W, leptons, 
IpT , and jets are incorporated. The effect of the underly- 
ing event, caused by interacting spectator quarks, is also 
included. The FMC program receives the particle four- 
momenta for each generated DYRAD event along with 
the next-to-leading order cross section prediction from 
DYRAD (which includes gluon loop effects) and produces 
kinematic distributions smeared by detector resolution, 
and sculpted by geometrical and kinematic acceptances 
and efficiencies. The FMC also reports event yield pre- 
dictions. The FMC successfully reproduces the kinematic 
features of inclusive W and Z boson production, as well 
as the features of vector boson production in association 
with a jet [8, 16]. 

For the W+jet data, we additionally require at least 



one "good" jet (£f * > 15 GeV and | < 2.4) that also 
passes the Aii]_j > 0.7 cut, where AR^ is the open- 



ing angle in ryiab-^iab space between the lepton and the 
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FIG. 2: The theoretical charged lepton <j> distribution in the 
Collins-Soper W rest-frame for the four regions, as gener- 
ated by DYRAD. The distributions are normalized to unity. 



TABLE II: The electron (N e ) and muon (JV M ) FMC event 
yields for inclusive W+jet production up to order a 2 . Con- 
tributions from backgrounds are not included. The uncer- 
tainties in the event yield predictions are dominated by the 
uncertainties associated with integrated luminosities, lepton 
identification and trigger efficiencies and the DYRAD predic- 
tion of the W+jet production cross section. Systematic un- 
certainties associated with PDF choice and Q 2 scale variation 
are not included. 



FMC event yields for inclusive W+jet production 


p¥ (GeV) 


N e 




15-25 
25-35 
35-65 
65-105 


3867 ± 137 
2632 ± 93 
2474 ± 87 
518 ± 18 


2027 ± 102 
1384 ± 66 
1314 ± 67 
279 ± 14 



leading "good" jet. The FMC event yields for inclusive 
IF+jet production up to order a 2 are presented in Table 
II. The Parton Distribution Functions (PDF) systemat- 
ics and the renormalization and factorization scale (Q 2 ) 
systematics will be included in Section VI. 

The FMC detector simulation, along with DYRAD, 
shows how the acceptances and efficiencies of the de- 
tector and the analysis cuts affect the </> distributions 
that are experimentally observed. Figure 3 shows the 
expected measurement of the <p distributions for the elec- 
tron dataset (the muon distributions are almost identi- 
cal) for the four pjf bins. The effects of the acceptances 
and efficiencies are significant; instead of two minima we 



FIG. 3: The expected CDF electron (j> distribution in the 
Collins-Soper W rest-frame for the four regions, after ex- 
perimental cuts and detector smearing, as generated by the 
FMC. The distributions are normalized to unity. The muon 
distributions are almost identical. 

observe two maxima. The main reason for this is the 
charged lepton and neutrino pr cuts, which limit the al- 
lowed (cos 9, 4*) phase space considerably. The FMC plots 
are normalized to the FMC signal event yields, and all 
experimental cuts have been applied. 



IV. ACCEPTANCES AND EFFICIENCIES 

The lepton identification and trigger efficiencies are 
measured by using the leptons from CDF Run la and 
lb Z data and by studying random cone distributions of 
leptonic W and Z decay Run la and lb data samples. The 
kinematic and geometrical acceptances are calculated us- 
ing the DYRAD event generator, which produces the SM 
prediction, and the FMC detector simulation, which pro- 
duces the CDF experimental expectation. 

We are especially interested in the product of overall 
acceptance times efficiency (ae) as a function of (cos 9,<j>) 
associated with each of the four pY^ bins We create 2- 
dimensional histograms of cos 9 vs. 4> f° r eacn of the four 
Pt bins, using the DYRAD simulation. This procedure 
is repeated after the events pass the FMC simulation, 
where the appropriate mixture of Run la and Run lb W 
leptons is used, based on FMC event yield predictions for 
all subdetectors. 

The resulting plots are shown in Figures 4, 5, and 6, 
for DYRAD, FMC-electrons, and FMC-muons respec- 
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FIG. 4: The cos 9 vs. <f> phase space for the four pt^ bins, for 
the DYRAD generator simulation (arbitrary units). 



FIG. 6: The cos 8 vs. <f> phase space for the four pt^ bins, for 
the muon FMC signal simulation (arbitrary units). 
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FIG. 5: The cos 6 vs. <j> phase space for the four pt^ bins, for 
the electron FMC signal simulation (arbitrary units). 



FIG. 7: Acceptance times efficiency for the electrons as a 
function of cos# and <f> in the Collins-Soper frame. 



tively. We subsequently divide the FMC 2-dimensional 
histograms by the corresponding DYRAD ones, produc- 
ing the 2-dimcnsional differential acceptance times effi- 
ciency ae(cos#, (f>) of Figures 7 and 8, for electrons and 



muons respectively. The overall acceptance times effi- 
ciency is higher for the electrons. These ae(cos9,cj)) val- 
ues are used for the cos 8- integration of the cross section, 
as described in Section VII. 
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FIG. 8: Acceptance times efficiency for the muons as a func- 
tion of cos d and <j> in the Collins-Soper frame. 



V. BACKGROUND ESTIMATION 

The main sources of background in the (W — ► ev e )+]ct 
and (W — > /if M )+jet processes are Z+jets events where 
the Z is misidentified as a W ("one-legged" Z), (W — > 
r^ T )+jet events, and QCD background resulting from jet 
misidentification. The W + 7 background is only a few 
events with run and event numbers the same as the W+"f 
events of [16]; it is treated as a systematic uncertainty, 
which also gives an indication of the radiative effects in 
the measurement. This uncertainty is very small (see 
Section IX). 

A small background contribution arises from ti pro- 
duction, where one of the produced W bosons decays 
lcptonically and the other W boson decays hadronically 
to jets. This background is estimated to be 30 ± 7 events 
for the electron sample [17] and 16 ± 3 events for the 
muon sample, a 0.3% effect. An equally small back- 
ground is the Z — > t + t~ production, where one of the 
tau leptons decays hadronically and the other one lepton- 
ically. This background is estimated to be 47 ± 1 events 
in the electron sample [17] , and 25± 1 events in the muon 
sample, a 0.5% effect. To demonstrate the insignificance 
of the ti and Z — > t + t~ backgrounds, wc perform our 
analysis including the charged lepton <j> distribution for 
these background events, in several possible shapes, for 
the four bins. The resulting change in the extracted 
values of the angular coefficients is negligible compared 
to our systematic and statistical uncertainties. Thus, we 
ignore the backgrounds associated with ti production and 



Z — ► t + t~ decays. 

Finally, the cosmic ray background in the muon IF+jet 
datasets is estimated to be significantly less than 0.1%, 
and is therefore neglected. 



A. One-legged Z background 

To study this background we generate a DYRAD sam- 
ple of Z+jet events and pass it through the FMC Monte 
Carlo simulation and the subsequent analysis program. 
This predicts how many Z bosons are misidentified as 
W bosons. In these cases, the Z bosons satisfy all 
kinematic and lepton identification cuts for W bosons, 
but one of their decay leptons, or legs, is undetected. 
The DYRAD cross section for iT+jet up to order a 2 is 
68.21 ± 0.37 pb. For this DYRAD simulation we used 
Q 2 = (Ml 010 ) 2 = 91.2 GeV, the CTEQ4M(A = 0.3 GeV) 
parton distribution functions, 0.7-cone jets, jet-jet angu- 
lar separation greater than 0.7 in r/i a b-0i a b space, and 
E 3 T > 10 GeV. At the FMC level, we impose our usual 
W boson event selection cuts and additionally require at 



least one "good" jet (E 3 ™ > 15 GeV and 



< 2.4) 



that also passes the AR]^ > 0.7 cut. These results 
are summarized in Table III. Overall we expect 123 ± 5 
electron one-legged-Z+jet events and 337 ± 18 muon one- 
legged-Z+jet events passing the W+jet cuts, without ap- 
plying any cut on the W transverse momentum. Compar- 
ing these numbers to the FMC event yields for VF+jet, 
the one-legged-Z+jet background is (1.14±0.06)% for the 
electron W+jet and (5.90 ± 0.43)% for the muon VF+jet 
sample. This background is higher for the muon sample, 
because of the limited coverage of the muon chambers, 
which is responsible for higher yields of one-legged muon 
Z bosons. 

To examine how this background affects the VF+jet 
lepton distribution, we plot the <p distribution for the 
leptons from these processes for the four bins (Fig- 
ures 9 and 10). We see that the same pattern of two max- 
ima at 77 and j s present. The background plots are 
normalized to the expected event yields from the FMC, 
multiplied by a factor of five (to make them visible) , and 
superimposed on the signal FMC distributions, normal- 
ized to the signal FMC event yields. We include the one- 
legged Z FMC <j) distribution in the complete theoretical 
prediction of the (f> distributions, in order to correctly 
extract the angular coefficients. 



B. (W — ► ri/)+jet background 

If the W boson decays to a t that subsequently decays 
lcptonically, the three final neutrinos contribute to the 
IfjT , which is incorrectly associated with a single neu- 
trino. The signal of one charged lepton along with the 
E T mimics that of a W directly decaying to the charged 
lepton. Most of the tau background is removed when we 
utilize the fact that the charged lepton and Et coming 
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FIG. 9: Electron cf> distributions for the four p T bins for 
electron W+jet FMC events (solid histogram) and for Z+jet 
FMC background (multiplied by 5), where one of the electrons 
from the Z decay is undetected and the other one passes the 
detection and analysis requirements (dashed histogram) . The 
histograms are normalized to the electron FMC signal event 
yields. 



FIG. 10: Muon <f> distributions for the four pif bins for muon 
W+jet FMC events (solid histogram) and for Z+jet FMC 
background (multiplied by 5), where one of the muons from 
the Z decay is undetected and the other one passes the de- 
tection and analysis requirements (dashed histogram). The 
histograms are normalized to the muon FMC signal event 
yields. 



TABLE III: Monte Carlo background estimation of the num- 
ber of electron and muon one-legged Z+jet events. The 
background fractions are calculated with respect to the FMC 
W+jet event yields. 



One-legged Z+jet background 


p™ (GeV) 




Fraction 




Fraction 


15-25 
25-35 
35-65 
65-105 


47 ± 2 
30 ± 1 
25 ± 1 
5 ± 


1.22 ± 0.07 % 
1.14 ± 0.05 % 
1.01 ± 0.05 % 
0.96 ± 0.03 % 


127 ± 7 
82 ± 4 
72 ± 4 
12 ± 1 


6.26 ± 0.47 % 
5.92 ± 0.40 % 
5.48 ± 0.41 % 
4.30 ± 0.42 % 



from the r decay are soft. As a result, the W transverse 
mass in the r events is significantly smaller than that in 
the electron or muon events. By applying the pr cuts for 
the leptons and the W transverse mass cut, we remove 
92% of the tau W+jet events at the DYRAD generator 
level. 

To study the remaining tau background we start 
with a tau W+jet DYRAD sample (Q 2 = (M^ 6 ) 2 , 
CTEQ4M(A = 0.3 GeV) parton distribution function 
and 0.7-cone jets in ryiatr^iab space), and we let the tau 
decay to an electron or a muon. We then vector-sum 
the three neutrinos resulting from the W and tau decays 
to form a single 1£t ■ Subsequently, we pass the events 
through the FMC detector simulator to see how many 



events pass the W+jet cuts after they are weighted by 
the detector acceptances and efficiencies. The branching 
ratios for the tau decays we use are 17.83 % for electrons 
and 17.37 % for muons [18]. At the FMC level, we require 
at least one "good" jet (E^ > 15 GeV and |r^* | < 2.4) 
that also passes the Ai?}^- > 0.7 cut. The tau back- 
ground results are presented in Table IV. Overall we 
expect 247 ± 9 tau electrons and 130 + 7 tau muons to 
infiltrate the W+jet samples, without applying any cut 
on the W transverse momentum. Comparing these num- 
bers to the FMC event yields for the electron and muon 
W+jet samples, the tau background is (2.28 + 0.12)% for 
the electron W+jet sample, and (2.28 ± 0.17)% for the 
muon W+jet sample. 

To see how this background affects the W+jet lepton 
<f> distribution, we plot the (j> distribution for the leptons 
resulting from leptonic tau decays in W+jets events for 
the four bins (Figures 11 and 12). We see that the 
pattern of two maxima at ^ and 4f is again present. The 
background plots are normalized to the expected event 
yields from the FMC, multiplied by a factor of five, and 
superimposed on the signal FMC distributions, normal- 
ized to the signal FMC event yields. We include the r- 
background FMC cf) distribution in the complete theoret- 
ical prediction of the </> distributions, in order to correctly 
extract the angular coefficients. 
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FIG. 11: Electron (f> distributions for the four bins for 
electron W+jet FMC events (solid histogram) and for tau 
W+jet FMC background (multiplied by 5), where the tau 
decays to an electron (dashed histogram). The histograms 
are normalized to the FMC signal event yields. 



FIG. 12: Muon <j) distributions for the four p!f. bins for muon 
W+jet FMC events (solid histogram) and for tau W+jet 
FMC background (multiplied by 5), where the tau decays to 
a muon (dashed histogram). The histograms are normalized 
to the FMC signal event yields. 



TABLE IV: Monte Carlo background estimation of the num- 
ber of electron and muon W+jet events, where the W decays 
to a tau and the electron or muon is the decay product of 
the tau. The fractions of the backgrounds are calculated with 
respect to the FMC W+jet event yields. 



(W — ♦ r v T — » v T v t fe/n e/^t)+jet background 


Pt (GeV) 




Fraction 




Fraction 


15-25 
25-35 
35-65 
65-105 


86 ± 3 
57 ± 2 
56 ± 2 
15 ± 1 


2.22 ± 0.11 % 
2.16 ± 0.10 % 
2.26 ± 0.11 % 
2.89 ± 0.22 % 


45 ± 2 
30 ± 2 
30 ± 2 
8 ± 


2.22 ± 0.15 % 
2.17 ± 0.18 % 
2.28 ± 0.19 % 
2.87 ± 0.14 % 



C. QCD background 

The QCD background in the case of inclusive W pro- 
duction and decay consists predominantly of dijet events, 
where one of the jets is misidentified as a lepton and the 
other one is not detected, resulting in the creation of 
$t ■ In the W+jet case, the QCD background is multi- 
jet events, where one of the jets is detected, one is lost or 
mismeasured (resulting in JpT ) and one is misidentified 
as a charged lepton to erroneously reconstruct a W. The 
number and distribution of QCD background events in 
the four pjf bins are determined from the Run la and 
Run lb CDF data. 

To measure the expected number of QCD background 
events in our data samples we look at leptons with isola- 



tion (ISO), defined in Section II, greater than 0.2. Our 
signal is in the ISO < 0.1 region and most of the events 
with lepton ISO > 0.2, but not all of them, are QCD 
background events. The upper histogram of Figure 13 
shows the isolation distribution of the electrons from 
W+jet events, for the first bin. When plotted on 
a semi-log scale, the ISO < 0.1 and the ISO > 0.2 re- 
gions can be approximated with two straight lines. The 
technique we use extrapolates the ISO > 0.2 line into the 
ISO < 0.1 signal region to calculate its integral and ob- 
tain the number of events in the signal region, using the 
assumption that the QCD background shape is not al- 
tered in that region. This method would give us the true 
number of QCD background events, if the ISO > 0.2 re- 
gion was filled exclusively with QCD events. In reality, 
only a fraction of these events are true QCD background, 
the rest being W+jet events. Since we expect to have 
some W+jet events in the region of lepton isolation from 
0.1 to 0.2, we fit the area above 0.2 with a straight line 
(in the semi-log histogram), which describes the QCD 
background. We also fit five continuous regions of lep- 
ton isolation, around the central region of ISO=0.20 to 
ISO=0.65 (namely 0.15-0.65, 0.25-0.65, 0.15-0.60, 0.20- 
0.65, and 0.25-0.70) to obtain a systematic uncertainty 
for this procedure. 

Since not all of the extrapolated region is QCD back- 
ground, we obtain a measurement of the percentage of 
the true QCD background in the electron W+jet sample 
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TABLE V: The linear least-squares fit parameterization of the A^l^ distribution for ISO > 0.2 electron W+jet events (second 
column) allows us to estimate the number of ISO > 0.2 W+jet events in the A^j!* > 2.5 region. The integral of this line, 
divided by the bin width (tt/30) of the A</>jt b j histogram, is the number of W+jet events with ISO > 0.2 and A^jt* > 2.5 
shown in the third column. These events are subtracted from the total number of (W+jet + QCD) background events in the 
ISO > 0.2 and Acfi^j > 2.5 region (fourth column). The result is divided by the total number of events in the ISO > 0.2 
region, to obtain an estimate of the fraction of true QCD background events (fifth column). 



p™ (GeV) 


Fit parameterization of electron 
W+jet events with ISO > 0.2 


Electron W+jet events with 
ISO > 0.2 and Ae^l^ > 2.5 


(W+jet)+QCD events with 
ISO > 0.2 and A^Jl^ > 2.5 


Fraction of 
true QCD events 


15-25 


2.21 x A^j^. - 0.05 


37.9 


257 


0.66=(257-37.9)/332 


25-35 


0.70 x &.4>Y-j + 0.78 


16.9 


98 


0.58=(98-16.9)/141 


35-65 


1.04 x A^* 13 .,. + 0.45 


20.6 


49 


0.31=(49-20.6)/91 


65-105 


0.17 x A0^ b . + 0.73 


7.5 


10 


0.13=(10-7.5)/20 



TABLE VI: The number of electron W+jet events, the number of QCD background events before correction and their percentage 
in the signal region, and the fraction of true QCD background events and their percentage in the signal region, for the four p% 
bins (see text for details). 



P? (GeV) 


Number of 
electron W+jet 


QCD events 
before correction 


Percentage of QCD 
before correction 


Fraction of 
true QCD events 


Percentage of 
QCD background 


15-25 
25-35 
35-65 
65-105 


5166 
3601 
3285 
624 


qco+26 
•J 00 -148 

14± 8 3 


8 - 18 i(V.81% 

9.80+']-I 4 % 
l-64l 4 ;' 6 7 °% 
2.24+^ 4 % 


66+ U U4 

58+ - 07 
u - oo -0.08 

1S+ 087 


5.40+^% 
5.68+° : «% 
0.51+J- 44 % 
0.29lJ-^% 




Electron Isolation (W+jet events, 15 < p y w< 25 GeV) 




Electron A<^ for ISO > 0.2 (W+jet events, 15 < p T W< 25 GeV) 



FIG. 13: Isolation of the electron in W+jet events (upper his- 
togram) and the difference in the azimuthal angles of the elec- 
tron and the jet for high isolation events (lower histogram). 
These histograms are used to estimate the QCD background 
in the electron W+jet data, for the 15 < p% < 25 GeV bin 
(see text). Corresponding histograms are used for the esti- 
mation of the QCD background in the three other electron 
W+jet p^ bins. 



above electron isolation of 0.1, by making a histogram of 
Acj}]^ for the events with ISO > 0.2, where is the 

difference in the <fi angle between the electron and the 
highest-i^T jet, with no other requirements for that jet. 
We expect the Acj)]^ distribution to be almost flat for 
the IF+jet events, because no correlation exists between 
the jet and the lepton <j) directions. In reality, this dis- 
tribution decreases at low Af^j^ , due to the application 
of the lepton isolation cut in our data. For QCD back- 
ground, we expect the A0]l^- between the highest Et 
jet and the jet resembling the lepton to peak at n. The 
lower histogram of Figure 13 shows the Acf^j for the 
events with lepton isolation greater than 0.2 for electron 
W+jet events and for the first bin. We fit the re- 
gion A^]^. < 2.5 (W+jet contribution) with a straight 

line. The region of the histogram A^]^ > 2.5 above 
that line corresponds to true QCD background. By di- 
viding this part of the histogram by the total number of 
events with ISO > 0.2, we determine the true fraction of 
QCD background in the ISO > 0.2 region. We expect 
the same fraction to be valid in the signal region (ISO 
< 0.1). Therefore, the number of true QCD background 
events is obtained by multiplying the number of ISO < 
0.1 events (as obtained by extrapolating the ISO > 0.2 
line into the signal region of the lepton isolation plot) by 
the QCD background fraction obtained from the Acj)]^ 
plot. The procedure is repeated for the four bins. Ta- 
ble V shows the extracted fraction of QCD background in 
the ISO > 0.2 region for the four pjf bins. The electron 
W+jet QCD background results are presented in Table 
VI. 



In the study of the QCD background in the muon sam- 
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Muon for ISO > 0.2 (W+jct events, 15 < p ? w< 25 GeV) 



FIG. 14: Isolation of the muon in W+jet events (upper his- 
togram) and the difference in the azimuthal angles of the 
muon and the jet for high isolation events (lower histogram). 
These histograms are used to estimate the QCD background 
in the muon W+jet data, for the 15 < < 25 GeV bin 
(see text). Corresponding histograms are used for the es- 
timation of the QCD background in the three other muon 
W+jet bins. At low A^J^ the distribution increases, due 
to bremsstrahlung processes associated with residual muon 
Z+jet background in the muon W+jet sample. These events 
disappear if we apply the zmuo_veto cut. (see text). 



pie we face a new problem. We originally apply a cut to 
muon W+jet data (zmuo_veto) in order to remove events 
that are consistent with the production of a Z boson, 
where one of the muons is non-isolated because it fails 
one (and only one) of the following cuts: 

• The muon isolation cut ISO < 0.1 

• The electromagnetic calorimeter cut £; EM < 2 GeV 

• The hadron calorimeter cut E RAD < 6 GeV 

These dimuon events are true Z bosons that look like W 
bosons because one muon does not pass one of the above 
cuts due to inner bremsstrahlung or bremmsstrahlung 
in the electromagnetic or hadronic calorimeters. The 
zmuo.veto cut mainly affects the tail of the muon iso- 
lation distribution (ISO > 0.2) and causes us to under- 
estimate the QCD background, since we use that tail to 
estimate it. Therefore, for the muon W+jet samples, for 
the purposes of determination of QCD background, we 
neglect this cut, in order to remove this bias at high muon 
isolation (ISO > 0.2) and make the transition from the 



low to high isolation smooth. Some of the muon Z back- 
ground is thus counted as QCD background; however we 
do not expect it to radically affect our QCD background 
estimation. In the isolation method we fit the back- 
ground starting from ISO=0.17 to ISO=0.40, to increase 
the statistical significance of our estimation. We also fit 
five continuous regions of lepton isolation, around the 
central region of ISO=0.17 to ISO=0.40 (namely 0.16- 
0.40, 0.18-0.40, 0.16-0.35, 0.17-0.40, and 0.18-0.45) to 
obtain a systematic uncertainty for this procedure. The 
upper histogram of Figure 14 shows the isolation distri- 
bution and fits for the muon W+jet events and for the 
first pf bin. 

We obtain a measurement of the percentage of the 
true QCD background in the muon W+jet sample above 
muon isolation of 0.2, by making a histogram of A<^]t_j 
for the events with ISO > 0.2, where A^j^ is the differ- 
ence in the <fr angle between the muon and the highest--&r 
jet, with no other requirements for that jet. The lower 
histogram of Figure 14 shows A</)|^ for the events with 
isolation greater than 0.2 for muon W+jet events, for 
the first pYf bin. The peak in the A(f>i-j = region is 
due to the muon bremsstrahlung processes that are not 
suppressed after we relax the zmuo.veto cut. We ignore 
these events when we fit to the straight line describing 
the W+jet events with high isolation muons. Table 
VII shows the extracted fraction of QCD background in 
the ISO > 0.2 region for the four pjf bins. The muon 
W+jet QCD background results are presented in Table 
VIII. For the highest muon p!^ bin the predicted number 
of true W+jet events is greater than the total number of 
events with ISO > 0.2 and A^} 3 ^- > 2.5, which results in 
a fraction of true QCD background events above ISO > 
0.2 equal to zero. 

After we calculate the percentage of the QCD back- 
ground in the signal region, we multiply it by the CDF 
W+jet event yields to obtain the absolute prediction of 
the number of QCD background events in each of the 
four p^ bins, for both and electron and muon W+jet 
data. The results are presented in Table IX. 

To complete the study of the QCD background we need 
to estimate its shape to properly include this background 
in the Standard Model prediction of the lepton <f> distri- 
bution in the CS frame, for each of the four pYj( bins. We 
plot cj) for the events with ISO > 0.2 and A^- > 2.5 
for the electrons and muon datasets, as shown in Figures 
15 and 16, respectively. We fit the distributions to the 
sum of two Gaussians and two straight lines. For the last 
p^ bin of the electrons and the last two p^ bins of the 
muons, there arc not enough statistics for the fit, so we 
use the total distributions (for 15 < pYp < 105 GeV) nor- 
malized to the number of events for those high p^ bins. 
We do not expect the shape of the QCD background to be 
significantly altered with increasing p^ . We assume that 
these distributions are the same as the ones in the signal 
region (ISO < 0.1) after they are properly normalized. 
We use these distributions to add the QCD background 
to the Standard Model prediction, after they are normal- 
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TABLE VII: The linear least-squares fit parameterization of the Ac^l^- distribution for ISO > 0.2 muon W+jet events (second 
column) allows us to estimate the number of ISO > 0.2 W+jet events in the A^j!* > 2.5 region. The integral of this line, 
divided by the bin width (tt/30) of the A^jt* histogram, is the number of W+jet events with ISO > 0.2 and A^jt* > 2.5 
shown in the third column. These events are subtracted from the total number of (W+jet + QCD) background events in the 
ISO > 0.2 and A0jLj > 2.5 region (fourth column). The result is divided by the total number of events in the ISO > 0.2 
region, to obtain an estimate of the fraction of true QCD background events (fifth column). 



Pt (GeV) 


Fit parameterization of muon 
W+jet events with ISO > 0.2 


Muon W+jet events with 
ISO > 0.2 and A^Jl^ > 2.5 


(W+jet)+QCD events with 
ISO > 0.2 and A^^ > 2.5 


Fraction of 
true QCD events 


15-25 


0.29 x Aci 1 ,^ + 1.31 


13 


164 


0.52=(164-13)/288 


25-35 


0.46 x Ar/>\ ab : j + 0.79 


12.8 


69 


0.40=(69-12.8)/140 


35-65 


0.25 x + 1.02 


10.5 


19 


0.14=(19-10.5)/61 


65-105 


x A^ + 1 


6.1 


1 


= 0/6 



TABLE VIII: The number of muon W+jet events without the application of the zmuo_veto cut, the number of QCD background 
events before correction and their percentage in the signal region, and the fraction of true QCD background events and their 
percentage in the signal region, for the four pf bins (see text for details) . 



P? (GeV) 


Number of 
muon W+jet 


QCD events 
before correction 


Percentage of QCD 
before correction 


Fraction of 
true QCD events 


Percentage of 
QCD background 


15-25 
25-35 
35-65 
65-105 


2779 
1943 
2002 
389 


280+^ 

103+i 29 

139 + |° 

11+° 
1J --o 


10.07+.^% 
5.30+°-^% 
6.94 + ^°% 
2.70+°;°°% 


52+ u - u5 

40+ - 08 

14+ - 22 

+1 
u -() 


5.24+^% 
2.12+1%% 
0.97+_lf r % 
0+f, 7 % 




Electron <|>, ISO > 0.2 and At^ > 2.5 (35 < p T W < 65 GeV) 




Electron t|), ISO > 0.2 and A$ H > 2.5 (65 < p T W < 105 GeV) 




Muon <>, ISO > 0.2 and A*, > 2.5 (15 < p T * < 25 GeV) 




Muon 4>, ISO > 0.2 and A$, > 2.5 (25 < p T W < 35 GeV) 




Muon t, ISO > 0.2 and A^> 2.5 (65 < p T W < 105 GeV) 



FIG. 15: The Collins-Soper <fi distribution of electrons from 
W+jet events with ISO > 0.2 and A^ 3 * > 2.5 for each of 
the four p% bins. These events are predominantly QCD back- 
ground events. We fit the distribution of the first three p% 
bins with two Gaussians on top of two straight lines. For 
the highest p% bin we use the distribution of the total QCD 
background, normalized to the number of the QCD events in 
this bin. (see text for details). 



FIG. 16: The Collins-Soper (j> distribution of muons from 
W+jet events with ISO > 0.2 and A^ > 2.5 for each 
of the four p^ bins. These events are predominantly QCD 
background events. We fit the distribution of the first two 
p% bins with two Gaussians on top of two straight lines. For 
the two highest p^ bins we use the distribution of the total 
background, normalized to the number of the QCD events in 
those bins. 
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TABLE IX: QCD background estimation for the electron and muon W+jet events. The fractions of the backgrounds are 
calculated with respect to the CDF Data W+jet events. 



QCD background 



Pt (GeV) 


N e 


Fraction 




Fraction 


15-25 
25-35 
35-65 
65-105 


279 tt L s 
205 tH 
17 ttl 

9 +12 
z -2 


5.40 t» «i % 
5.68 +0 2 H % 
0.51 +£■** % 
0.29 +IH % 


148 
40 +g° 
18 tT 8 

o tl° 


5.24 ±J;£ % 
2.12 t 2 ;S % 
0.97 tlf 7 % 

o +g- 7 % 



TABLE X: Summary of electron W+jet backgrounds. The background fractions are calculated with respect to the FMC signal 
event yields for the electroweak backgrounds and with respect to the data for the QCD background. 



Electron W+jet Backgrounds 


Background 


Pt' =15-25 GeV 


Pt' =25-35 GeV 


p? 7 =35-65 GeV 


p^ =65-105 GeV 


W -> TV T 

Z -> e+e" 
QCD 


86 + 3 (2.22 %) 
47 + 2 (1.22 %) 
279+H (5.40 %) 


57 + 2 (2.16 %) 
30 ± 1 (1.14 %) 
205+.^ (5.68 %) 


56 + 2 (2.26 %) 
25 ± 1 (1.01 %) 
17+tI (0.51 %) 


15 ± 1 (2.89%) 
5 + (0.96 %) 
2+ 2 2 (0.29 %) 



TABLE XI: Summary of muon W+jet backgrounds. The background fractions are calculated with respect to the FMC signal 
event yields for the electroweak backgrounds and with respect to the data for the QCD background. 



Muon W+jet Backgrounds 


Background 


Pt 1 =15-25 GeV 


Pt~ =25-35 GeV 


=35-65 GeV 


Pt* =65-105 GeV 


W -> TV T 

Z -> v+pT 
QCD 


45 + 2 (2.22 %) 
127 ± 7 (6.26 %) 
148+^ (5.24 %) 


30 + 2 (2.17 %) 
82 + 4 (5.92 %) 
40+f (2.12 %) 


30 + 2 (2.28 %) 
72 + 4 (5.48 %) 
18+^ (0.97 %) 


8 + (2.87%) 
12 ± 1 (4.30 %) 
0+o° (0 %) 



TABLE XII: The expected total event yields for inclusive W+jet production. The signal and electroweak backgrounds are 
calculated up to order a 2 . The PDF and Q 2 systematics have also been included (second set of uncertainties). 



Expected signal+background event yields for inclusive W+jet production 


(GeV) 


Electrons 


Muons 


N e (Signal) 


N e (Background) 


N e (Total prediction) 


N„ (Signal) 


JV M (Background) 


JV M (Total prediction) 


15-25 
25-35 
35-65 
65-105 


3867 ± 137 
2632 ± 93 
2474 ± 87 
518 ± 18 


292tf 
98+f[ 
22+1 2 


42791^ +880/-660 
2924li°°+598/-408 
2572+i° 2 +562/-383 
54of 22 ,+118/-81 


2027 ± 102 
1384 ± 66 
1314 ± 67 
279 ± 14 


320+^ 
152+^0 

120l?l 

20t\° 


2347t^ 4 +484/-330 
1536+^2 +329/-224 
1434+_™+312/-213 

2 "ii5+ 66 /- 45 



ized to the expected number of QCD background events, 
given by Table IX. 



D. Summary of backgrounds and Standard Model 
event yields prediction. 

Backgrounds for electron and muon W+jet events for 
each of the four pYp bins arc summarized in Tables X 
and XI respectively. We obtain the total W+jet event 
yield prediction by adding these backgrounds to the FMC 
W+jet signal prediction of Table II. To obtain the final 
uncertainties, we add linearly the uncertainties associ- 
ated with the W+jet signal and electroweak background 
and add the result to the QCD background uncertainty 
in quadrature. The total W+jet event yields after the 
inclusion of the backgrounds are presented in Table XII. 
The PDF and Q 2 systematic uncertainties are also in- 



cluded (see Section IX). 

VI. COMPARISON BETWEEN EXPECTED 
AND OBSERVED W DISTRIBUTIONS 

We study the expected (FMC) W kinematical distri- 
butions after the inclusion of backgrounds and compare 
them to the experimental distributions. Figures 17 and 
19 show the W transverse momentum for electrons and 
muons respectively. The observed and simulated distri- 
butions have been normalized to unity. We observe good 
agreement between the observed and simulated pYp distri- 
butions. Figure 18 shows the W transverse mass distribu- 
tion for the electron W+jet datasct and for the DYRAD 
events passed through the FMC detector simulation. Fig- 
ure 20 shows the same distributions for the four bins. 
Figures 21 and 22 show the same distributions for the 
muon W+jet datasets. The observed and simulated dis- 
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20 40 60 80 100 20 40 60 80 100 120 

p T * (GeV) (Electron) M T W (GeV) (Electron, p T w > 1 5 GeV) 



FIG. 17: The transverse momentum of the W for the elec- 
tron W+jet data sample (points) along with the FMC sig- 
nal simulation including backgrounds (solid histogram). The 
backgrounds are multiplied by 5, to be visible. The data and 
expected signal+background distributions are normalized to 
unity. 



FIG. 19: The transverse momentum of the W for the muon 
W+jet data sample (points) along with the FMC signal sim- 
ulation including backgrounds (solid histogram). The back- 
grounds are multiplied by 5, to be visible. The data and 
expected signal+background distributions are normalized to 
unity. 
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FIG. 18: The transverse mass of the W for the electron W+jet 
data sample (points) along with the FMC signal simulation 
including backgrounds (solid histogram). The backgrounds 
are multiplied by 5, to be visible. The data and expected 
signal+background distributions are normalized to unity. 



FIG. 20: The transverse mass of the W for the electron W+jet 
data sample (points) along with the FMC signal simulation 
including backgrounds (histogram), for the four bins. The 
data and expected signal+background distributions are nor- 
malized to unity. 
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(Muon, 25 < p T w < 35 GeV) 




(Muon, 35 < p T w < 65 GeV) 
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FIG. 21: The transverse mass of the W for the muon IF+jet 
data sample (points) along with the FMC signal simulation 
including backgrounds (solid histogram). The backgrounds 
are multiplied by 5, to be visible. The data and expected 
signal+background distributions are normalized to unity. 



FIG. 22: The transverse mass of the W for the muon IF+jet 
data sample (points) along with the FMC signal simulation 
including backgrounds (histogram), for the four bins. The 
data and expected signal+background distributions are nor- 
malized to unity. 



tributions are again normalized to unity. In all of the 
above plots, the FMC distributions are produced with 
properly weighted signal and background contributions, 
for electron and muon detector regions. 

VII. DIRECT MEASUREMENT OF THE 
AZIMUTHAL ANGLE OF THE CHARGED 
LEPTONS FROM W DECAYS IN THE 
COLLINS-SOPER FRAME 

For each W event we boost to the W rest-frame to 
calculate the azimuthal angle of the charged lepton. The 
longitudinal momentum of the W {p^) is not known, 
because the longitudinal momentum of the neutrino is 
not measurable, so we use the mass of the W to constrain 
it. For a particular event, the longitudinal momentum 
of the neutrino is constrained by the mass of the W, 
according to the equation: 

P ' = fa ± El \/ A2 - A ^ 2 (Pt) 2 ) - ( 2 ) 

where 

A = M 2 W + (p^y - [ p i T f - {f T )\ (3) 

E l is the energy of the charged lepton, p l T is its transverse 
momentum, p l z is its longitudinal momentum, pi^ is the 
neutrino transverse momentum, and p™ is the transverse 



momentum of the W . This equation is unique for every 
event, since the kinematics of the lepton and neutrino, as 
well as the mass of the W, contribute to the shape of the 
curve p v z = f(Mw)- If the mass of the W was known on 
an event by event basis, there would be a two-fold ambi- 
guity in the value of p v z of the neutrino in the laboratory 
frame. Because the W boson has a finite width given by 
a PDF-convoluted Breit-Wigner distribution, BW(M), 
we actually have two distributions of possible values of 
p" z , BW{M{p v z )), where M{p u z ) is the mass of the W as a 
function of the neutrino longitudinal momentum for the 
particular kinematics of the event. 

The choice of one of the two neutrino longitudinal mo- 
mentum solutions does not affect the <j) analysis, since 
both solutions result in the same charged lepton <ft in the 
CS frame. For this analysis, only the choice of the W 
mass is of interest. The choice is made based on the 2- 
dimensional M w vs. Mjf histo grams constructed with 
DYRAD events. For a specific we use a probabil- 
ity distribution of W masses and randomly select one for 
each event, based on that distribution. This method was 
devised to better reconstruct the | cos((9)| distribution [8], 
since the polar angle is very sensitive to the selection of 
the W mass. In our analysis, the azimuthal angle is not 
affected by the choice of mass, so the answer is almost 
the same even if we choose a mass based on the Brcit- 
Wigner distribution and the requirement that the mass 
is greater than the measured transverse mass. 

After obtaining a <j) for every event, we proceed to an- 
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alyze our sample. Theoretically, the W differential cross 
section, integrated over cos 9 and y is given by: 



where 



^°" 2 j2 = C(l + Pi cos <j> + f3 2 cos 2<j) 



d{pW) 2 d<j) 



+/3 3 sin0 + /3 4 sin2</>), 



(4) 



where 



C 



da 



2tt d(p^ 



W\2 ' 



01 



^ Ar(p¥) ,^^fl. (5 ) 



The theoretical <j> distributions for the charged lepton 
from W boson decay in W+jet production are shown in 
Figure 2. 

From Equations (4) and (5), the reader might con- 
clude that only the A 2 , A3, A$, and A-j coefficients are 
measurable with the cj> analysis, since the other angular 
coefficients are integrated out. However, in the actual 
W+jet data samples, what we measure is the number of 
events: 



N(p$,<t>) = 



da 



d{p^) 2 d<j)d cos 9 



ae(pr , cos 9, (f>)dcos9 J Ldt 



N bg {p w 



T ) 1 



(6) 



where C is the instantaneous luminosity and 
ae(pj? , cos 9, 4>) is the overall acceptance times effi- 
ciency, determined in Section IV, for a particular W 
transverse momentum and region in the (cos9,(j)) phase 
space. The quantity Nt, g (p^ , 4>) is the background for 
the given <j> bin and p? , estimated in Section V. Com- 
bining Equations (6) and (1), the measured distribution 
is 



N{p™,4>) = C'(f- 1 (p™,<P)+J2MPT)fi(PT,<P)) 



i=0 



w 



+ N bg (p™ 



(7) 



where C — C J Ldt. The /j are fitting functions, which 
are integrals of the product of the explicit functions 
gi(cos9,4>) and ae(cos9,<p): 



/ gi(0,(j>)ae(pT , cos 9, <j))d cos 9, (8) 
Jo 



0-1 (M 


) = l + cos 2 <9 


.90 (M 


) = ~(1 - 3cos 2 6») 
2 


9i(0,(p 


) = sin 29 cos 


.92 (9,cf> 


) = i sin 2 9 cos 2(f> 


93(0,0 


) = sin 9 cos 


9^(0,0 


) = cos 9 


95(0 A 


) = sin 2 9 sin 2<j) 


96(0,<j> 


) = sin 29 sin 


97(0,0 


) = sin 9 sin 



(9) 



Because 



cos( 



-1.....7 



we multiply the gi(9,0) functions by 
, 0) before integrating over cos 9, no /j is ex- 
actly zero and all of the angular coefficients Ai are in 
principle measurable. We have verified that the FMC- 
simulated distributions, fitted with a linear combina- 
tion of result in angular coefficient values consistent 
with the SM predictions [7]. This result supports the 
self-consistency of the method. 

We use Simpson integration for the calculation of the 
/, fitting functions given by Equation (8). The explicit 
functions 9i(9, 0) are integrated over cos#, after they are 
weighted with the value of ae(p^ , cos 9, 0) extracted from 
the 2-dimensional histograms of Figures 7 and 8. 

Although the use of Equation (7) allows us in princi- 
ple to measure all of the angular coefficients, in reality, 
the current statistics do not allow us to make a signifi- 
cant measurement of angular coefficients other than A 2 
and A 3 . This is due to the fact that the fitting functions 
/i#2.3 are small, and the distributions are insensitive 
to large variations of the corresponding angular coeffi- 
cients. Figure 23 shows how the expected electron dis- 
tributions are modified as the angular coefficients Ai are 
varied, one coefficient at a time (the muon distributions 
are almost identical). Using Equation (7), we vary Ao, 
A 2 , and A 3 from to 1 with a step size of 0.1, and A4 
from to 2 with a step size of 0.2. We find that only A 2 
and A3 strongly affect the azimuthal distributions, thus 
only these two angular coefficients are measurable with 
the our current analysis. Large variations of A and A 4 
result in small changes in the distributions, hence the 
uncertainties associated with the measurement of Aq and 
A 4 are large; these two coefficients cannot be measured 
in a statistically significant manner with the analysis. 
The same is true for A\, A 5 , Aq, and A-j, all of which are 
consistent with zero for our current experimental preci- 
sion. 

Figure 24 shows the observed CS electron distribu- 
tions for CDF electron VF+jet data for the four pjf bins. 
Figure 25 shows the corresponding distributions of the 
CDF muon VF+jct data. The solid lines are the SM the- 
oretical predictions including backgrounds, whereas the 
points correspond to CDF VF+jet data (the error bars 
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FIG. 23: The predicted electron (j> distributions in the Collins- 
Soper W rest-frame for the first pr bin, varying only one 
angular coefficient at a time and keeping the other angular 
coefficients at their Standard Model values. Only Ai and 
Az significantly affect the shape of these distributions. The 
same is true for the other p% bins and also for the muon <f> 
distributions. 



FIG. 24: The <p distributions for the electron CDF data 
(points) the SM Monte Carlo (solid lines) and the result of 
the fit (dashed lines) for the four bins of p% . The errors are 
only statistical. The fit is performed from 7r/2 to 37r/2 and 
resulted in x 2 / dof equal to 2.32, 1.80, 2.18, and f .1 1 for the 
four bins respectively (11 degrees of freedom). All distribu- 
tions are normalized to unity. 



are statistical only) . The theoretical prediction for the <j> 
distributions is constructed using Equations (7) and (8). 
The free parameters are the angular coefficients Ai. The 
background <f> shapes are given by Figures 9, 11, and 15, 
for electrons and Figures 10, 12, and 16 for muons, nor- 
malized to the event yields of Tables X and XI respec- 
tively. The expected signal is normalized to the FMC 
signal event yields of Table XII, and subsequently the 
backgrounds are added to construct Nip™ ,<j>), accord- 
ing to Equation (7). The total theoretically predicted 
distributions along with the experimental ones, are fi- 
nally normalized to unity. The experimental results are 
in good agreement with the Standard Model prediction, 
which includes the effects of W polarization and QCD 
contributions up to order a 2 a . 



VIII. MEASUREMENT OF THE ANGULAR 
COEFFICIENTS 

The values of the angular coefficients A 2 and A 3 are 
extracted using the least-squares fitting method and the 
data associated with Figures 24 and 25. The least- 
squares fit is performed over the negative x-axis of the CS 
frame (tt/2 < <p < Sir/ 2) for the following two reasons. 

Firstly, if a single jet perfectly balances the W boson, 



its momentum will be placed on the positive x-axis in 
the CS frame. In reality, the leading jet will be in the 
x > region of the z — x plane, in proximity to the £-axis, 
as seen in Figure 26 for the electron W+jet data. The 
leading jet's <j) in the CS frame will almost always be less 
than 7r/2 or greater than 37r/2, as shown in Figure 26. A 
kinematic correlation exists between the angular separa- 
tion Ai? between the jet and the lepton in the </>i a b — /?i a b 
space and the CS cj> of the lepton, as shown in Figure 26. 
The situation is similar for the possible subleading jets in 
the VF+jet events (Figure 27). VF+jet events with more 
than two jets are not modeled in DYRAD simulation; 
their presence in the data creates extra biases in the low 
and high regions of the lepton <f> distributions. Because 
of the lepton-jet angular separation and lepton isolation 
requirements in our T'F+jet datasets we obtain a bias- free 
measurement of the angular coefficients A 2 and A 3 if we 
exclude the positive- x half-plane region of the CS frame. 

Secondly, the term A 3 f 3 {cf)) in Equation 7 is the small- 
est measurable term with our data. Therefore, a more 
significant measurement of the angular coefficient A 3 is 
obtained in the CS <f> region where the rest of the terms 
(and mainly the predominant A^f^if) term), contribute 
less. The ratio A 3 f 3 ((p) / A±f±{<f>) is significantly larger 
in the 7r/2 < <f> < 3tt/2 region, and thus a more sensi- 
tive measurement of A 3 is obtained in this region. We 
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FIG. 25: The distributions for the muon CDF data (points) 
the SM Monte Carlo (solid lines) and the result of the fit 
(dashed lines) for the four bins of p^ . The errors are only 
statistical. The fit is performed from 7r/2 to 37r/2 and resulted 
in x 2 /dof equal to 0.41, 1.39, 1.33, and 1.71 for the four bins 
respectively (11 degrees of freedom). All distributions are 
normalized to unity. 



FIG. 26: The p z vs. p x , p y vs. p x , and cj> of the leading jet 
in the CS frame and the AR between the jet and the lepton 
in the laboratory frame vs. lepton <j> in the CS frame, for 
electron data. 



normalize the theory to data from 7r/2 to 3tt/2 before we 
start the fitting procedure, which is carried out in the 
x < region of the z — x plane. 

We use the MINUIT % 2 minimization program [19] to 
fit the electron and muon (f> distributions to the fitting 
functions fi. Since these functions are not linearly in- 
dependent, we cannot fit with all parameters free. For 
this reason we keep the angular coefficients Aq and 
fixed at their SM values and allow A 2 and A 3 to vary. 
After we extract values for A 2 and ^3, we fix these coef- 
ficients at these values, and we repeat the fit procedure 
varying only the Aq and A4 angular coefficients. The 
angular coefficients A\, A5, Aq, and Aj are always kept 
fixed at their SM values, since the theoretical prediction 
for these coefficients is very close to zero and the varia- 
tion for the first 100 GeV of pj? is small in comparison 
to the experimental precision. We expect large statistical 
uncertainties for the extracted values of A and A4, since 
they do not significantly affect the <f> distribution. Large 
variations in their value only slightly alter the leptons' (j) 
angular distribution. 

The results of the MINUIT fits are shown as dashed 
histograms in Figure 24 for the electron W+jet data 
and Figure 25 for the muon VF+jet data. Our mea- 
surements of the angular coefficients for the electron and 
muon IF+jet data are presented in Figures 28 and 29 



GeV GeV 




FIG. 27: The p z vs. p x , p y vs. p x , and <f> of the extra jets 
in the CS frame and the AT? between the jets and the lepton 
in the laboratory frame vs. lepton (j> in the CS frame, for 
electron data. 
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FIG. 28: The measurement of the angular coefficients for the 
W+jet electron data (points) and the SM prediction up to 
order o? s (line). The errors are only statistical. 



FIG. 29: The measurement of the angular coefficients for the 
VF+jet muon data (points) and the SM prediction up to order 
«j (line). The errors are only statistical. 



respectively. The bin centers are determined using the 
average value of pjf for the range of the four pYp bins. 
The measured angular coefficients associated with the 
electron and muon VF+jet data agree with the SM pre- 
diction and with each other. We emphasize that the SM 
prediction is only up to order o? a in QCD. The statistical 
uncertainties for Aq and A4 are very large, as expected, 
making the measurement of these coefficients unrealistic 
using the azimuthal angle analysis. 

Assuming weak-interaction lcpton-universality, we 
combine the measurements of the angular coefficients 
obtained from the electron and muon H^+jet datasets, 
treating them as the results of two separate experi- 
ments. If A e and A^ are the electron and muon mea- 
surements with statistical uncertainties u e and re- 
spectively, then the combined measurement is A com b = 
(Ae/er? + ^/o-^)/(l/(7g + 1/tr^), with statistical uncer- 
tainty (7 com b = (l/oi + l/er^) -1 / 2 . The result of this 
statistical combination, along with the SM prediction, is 
presented in Figure 30. 



IX. SYSTEMATIC UNCERTAINTIES 
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The systematic uncertainties associated with the mea- 
surement of the angular coefficients Ai are related to the 
jet definition and energy scale, the selection of the W 
mass on an event-by-event basis, the background estima- 
tion, possible presence of W + 7 events in our datasets, 



FIG. 30: The measurement of the angular coefficients for the 
combination of electrons and muons (points) and the SM pre- 
diction up to order o? a (line). The errors are only statistical. 



23 



the assumed values of Aq and A4, the choice of parton 
distribution functions, and the renormalization and fac- 
torization scale Q 2 of the event. The jet systematic un- 
certainties, the variation of the A and A4 values, and Q 2 
scale uncertainty are the dominant sources of systemat- 
ics. 



A. Jet systematic uncertainties 

The number of data events passing the jet cuts is af- 
fected by the systematic uncertainties associated with the 
jet Et scale and the rapidity requirement. The same sys- 
tematic uncertainty has an effect on the measurement of 
the angular coefficients. 

The uncertainty on jet Et scale depends on the 
calorimeter stability, relative energy scale corrections, ex- 
tra interactions, and underlying event corrections. The 
total uncertainty is a quadratic sum of these effects. The 
systematic uncertainty in the jet energy scale affects the 
reconstruction of the IpT and the W boson. For ev- 
ery FMC VF+jet event, we shift the energy of the jet 
by a _|_ = 85%VE&, where E' et is the energy of the jet 
in GeV, without changing its direction. We then cor- 
rect the $t value and recalculate all the kinematic vari- 
ables associated with the W boson, jet, and $t ■ We 
subsequently extract the new acceptance times efficiency 
ae(cos 0, </>) and analyze the data. We repeat this pro- 
cedure for the energy shifted by er_ = — 85%V E^ ct and 
calculate the systematic effect of the jet energy scale on 
the measurement of the angular coefficients, presented in 
Table XIV for the electron, muon, and the combination of 
the two results. To obtain the combined results, wc com- 
bine the electron and muon measurements for each p^ 
bin and for each choice of energy shift, using the sta- 
tistical uncertainties of the central measurements. The 
difference between the shifted combined values and the 
central combined value determines the systematic uncer- 
tainty on the combined measurement. The same method 
is used for all the systematic uncertainty estimates. 

We vary the jet E T cut by {6E T ) = ±850 MeV in both 
data and MC and repeat the analysis each time, to deter- 
mine its effect on the measurement of the angular coeffi- 
cients and on the FMC-prediction of signal event yields. 
Table XIII shows the systematic uncertainty in the mea- 
surement of the W+jet event yields associated with the 
jet Et cut variation, for the four p^ bins. Overall, there 
is a +6.4%/ — 5.8% effect in the electron event yields and 
a +6.0%/ — 5.7% effect in the muon event yields due to 
the jet Et cut. 

The uncertainty on the rapidity r/i^ of the jet is Srji^b = 
±0.2. We vary the jet rjub cut from 2.2 to 2.6 to obtain 
the variation in the data event yields presented in Table 
XIII, for the four p™ bins. Overall, there is a +2.2%/ - 
2.4% effect in the electron event yields and a +0.7%/ — 
2.2% effect in the muon event yields due to the jet 77^ 
cut. 



In order to obtain an estimate of the systematic un- 
certainty in the measurement of the angular coefficients 
associated with the jet Et and ?7i a b cuts, we run the anal- 
ysis for 11 values of the E^ cut, from 14.15 GeV to 15.85 
GeV, and for five values of \r^ h \ cut; from 2.2 to 2.6. We 
record the variations in the measurement of the angular 
coefficients for electrons, muons, and the combination of 
the two results. The results for the four pYf bins arc 
presented in Table XIV. 



B. Systematic uncertainty due to W mass selection 

As previously discussed, in order to boost to the W 
rest-frame, a mass value is selected for the W boson. We 
have four different methods for selecting this mass on an 
event-by-event basis. We investigate how each mass se- 
lection method affects our angular coefficients measure- 
ment. The first method selects a Brcit-Wigner mass, 
which is greater than the measured transverse mass of 
the W boson. The second method selects the greater 
of the pole mass and the transverse mass. In the third 
method we select the pole mass or, in case it is less 
than the transverse mass, we select a Brcit-Wigner mass, 
which is greater than the transverse mass. Finally the 
fourth method (default) selects a mass based on the dis- 
tribution that results from the slice of the theoretical 
(DYRAD) M w vs. M^f 2-dimensional histogram (for 
IF+jet events) at the measured transverse mass of the 
W boson. This last method is preferred because it re- 
moves some biases in the measurement of the polar angle 
0. In the <j> analysis, the systematic uncertainty on the 
azimuthal angle cf) due to the selection of the mass of the 
W is minimal. We run the analysis for the four mass 
selection methods and record the variations in the mea- 
surement of the angular coefficients for electrons, muons, 
and the combination of the two results. All methods give 
almost identical measurements of <f>. The systematic un- 
certainties for the four p^ bins are presented in Table 
XIV. 



C. Backgrounds estimate systematic uncertainty 

There is an uncertainty in the estimation of the back- 
grounds, given by the uncertainties in Tables X and XI. 
We vary our prediction from the highest value to the 
lowest possible value for every background as well as the 
FMC signal event yields. These uncertainties do not in- 
clude the PDF and Q 2 systematics. For each variation, 
we rerun the analysis programs for the electron and muon 
case, and we also combine the results. The systematic 
uncertainties are presented for the four p^ bins in Table 
XIV. 
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TABLE XIII: Systematic uncertainties on VK+jet event yields due to the t and f^ b cuts, and the total systematic uncertainty 
due to these sources, for the four bins. 



VF+jet event yield systematic uncertainties due to the E^f ' and 77^^ cuts 


Pt (GeV) 


Charged Lepton 


SN due to 


SN due to rjlH 


Total systematic uncertainty 


15-25 


electron 
muon 


-515 
+308 
-290 


+ lbb 

-154 
+34 

-81 


+613 

-538 
+310 

-301 


25-35 


electron 
muon 


+i7y 

-175 
+95 
-79 


+78 

-104 
+12 
-44 


+195 

-204 

+y6 

-90 


35-65 


electron 
muon 


+38 

— 50 
+ 12 
-29 


+bU 

-43 

-29 


+63 

-66 
+ 12 
-41 


65-105 


electron 
muon 


+0 

-1 

+0 

-0 


+1 

-2 
+0 
-0 


+ 1 

-2 
+0 
-0 



D. W + 7 systematic uncertainty 

The IF+jet angular distribution can be affected by 
TF+7 production, for a hard 7 well-separated from the 
charged lepton from the W decay. Some of the events in 
our datasets are consistent with W + 7 production, ac- 
cording to [16]. We remove those events and remeasure 
A2 and ^3. The variation from the original measurement 
is treated as a systematic uncertainty. The systematic 
uncertainties for the four pj? bins are presented in Table 
XIV. 



E. A and A4 variation systematic uncertainty 

In our analysis we keep A and A 2 fixed at their SM 
values. To check how this affects our measurement, we 
set Aq and A 2 at minimum and maximum values in all 
possible combinations and repeat the analysis four times 
(Ao(min)=0, ^4o( m ax)=l, A4(min)=0 and A4(max)=2). 
The systematic uncertainties for the four pjf bins are 
presented in Table XIV. 



transverse momentum of the W , instead of the default 
square of the pole mass of the W boson. The systematic 
uncertainties for the four p^ bins are presented in Ta- 
ble XIV. If we try all Q 2 choices provided by DYRAD 
(total invariant mass squared, dynamic mass squared, 
total energy of the W squared, and transverse energy 
of the leading jet squared, in addition to the two men- 
tioned above), we end up with a systematic uncertainty 
of +19% / -10% on the DYRAD cross section, which 
affects both the central FMC signal event yields and the 
electroweak backgrounds. These variations are used for 
the estimation of the total FMC event yields systematic 
uncertainty due to Q 2 scale variation. 

H. Overall analysis systematic uncertainties 

Table XV summarizes the total systematic uncertain- 
ties for the A 2 and A 3 measurement, for the four p)^ 
bins and for the electron, muon, and combined results. 
To populate this table, we combine the systematics de- 
scribed above and presented in Table XIV. 



F. PDF systematic uncertainty 

To study the uncertainty associated with the parton 
distribution functions, we use the MRS A' \a s (Mz) — 
0.105 and A = 0.150] PDF [20] and repeat the analysis. 
The systematic uncertainties for the four p^ bins arc 
presented in Table XIV. When we use all PDFs of the 
MRSA and CTEQ families, we end up with a systematic 
uncertainty of ±11% on the DYRAD cross section, which 
affects both the central FMC signal event yields and the 
electroweak backgrounds. These variations are used for 
the estimation of the total FMC event yields systematic 
uncertainty due to choice of PDF. 

G. Q 2 Systematic uncertainty 

Finally we change the renormalization and factoriza- 
tion scale Q 2 so that it is equal to the square of the 



I. Overall systematic uncertainties in data and 
Monte Carlo event yields 

Combining the data event yield systematics due to 
Ej? and rfZu cut variations in quadrature, we get the 
final data event yields presented in Table XVI. Com- 
paring with the FMC event yields of Table XII, we see 
that there is a reasonable agreement with the SM predic- 
tion. In Table XII we have also included the PDF and 
Q 2 FMC systematic uncertainties described earlier, com- 
bined in quadrature to give a systematic uncertainty of 
+22% / —15% on the FMC signal event yields and elec- 
troweak background. We do not expect perfect agree- 
ment since the DYRAD generator produces up to two 
jets with E^ > 10 GeV (order a 2 ), while in the data we 
have many events with more than two jets with E^ * > 10 
GeV. If we impose a cut on the number of jets in the CDF 
data, by not accepting more than two jets in an event, 
and applying strict cuts on at least one jet, the disagree- 



25 



TABLE XIV: Systematic uncertainties in the measurement of A2 and A3 along with their sources, for electron and muon W+jet 
events and the combination of the electron and muon results. 



Source of systematic uncertainty 


p¥ (GeV) 


Electrons 


Muons 


Combination 


8A a 


SA 3 


SA 2 


6A3 


SA 2 


SA 3 


Jet Et cut 


15-25 
25-35 
35-65 
65-105 


+U.IWY6 

-0.0533 
+0.0395 
-0.0109 
+0.0433 
-0.0000 
+0.0000 
-0.0617 


+0.0000 

-0.0064 
+0.0047 
-0.0039 
+0.0030 
-0.0033 
+0.0167 
-0.0153 


+0.0146 

-0.0502 
+0.0676 
-0.0559 
+0.0177 
-0.0131 
+0.1519 
-0.0000 


+0.0036 

-0.0032 
+0.0073 
-0.0017 
+0.0042 
-0.0011 
+0.0055 
-0.0808 


+0.0230 

-0.0523 
+0.0490 
-0.0261 
+0.0357 
-0.0000 
+0.0556 
-0.0299 


+0.0000 

-0.0043 
+0.0050 
-0.0027 
+0.0031 
-0.0022 
+0.0024 
-0.(3374 


Jet ryiab cut 


15-25 
25-35 
35-65 
65-105 


+0.0043 

-0.0324 
+0.0007 
-0.0346 
+0.0093 
-0.0049 
+0.0000 
-0.0617 


+0.0000 

-0.0036 
+0.0012 
-0.0022 
+0.0023 
-0.0005 
+0.0022 
-0.0052 


+0.0098 

-0.0374 
+0.0299 
-0.0221 
+0.0070 
-0.0075 
+0.1021 
-0.0000 


+0.0044 

-0.0000 
+0.0000 
-0.0085 
+0.0038 
-0.0028 
+0.0000 
-0.0557 


+0.0002 

-0.0175 
+0.0106 
-0.0304 
+0.0070 
-0.0032 
+0.0112 
-0.0258 


+0.0000 

-0.(3018 
+0.00(30 
-0.0032 
+0.0013 
-0.0000 
+0.0000 
-0.(3203 


Jet energy scale 


15-25 
25-35 
35-65 
65-105 


+0.0594 

-0.0588 
+0.0851 
-0.1105 
+0.0017 
-0.0514 
+0.2468 
-0.3821 


+0.0042 

-0.0014 
+0.0096 
-0.0020 
+0.0204 
-0.0069 
+0.0836 
-0.0489 


+0.0499 

-0.0694 
+0.0860 
-0.0994 
+0.0142 
-0.0850 
+0.2989 
-0.3186 


+0.0075 

-0.0044 
+0.0261 
-0.0141 
+0.0183 
-0.0102 
+0.0382 
-0.0176 


+0.0560 

-0.0626 
+0.0854 
-0.1068 
+0.0054 
-0.0614 
+0.2660 
-0.3587 


+0.0052 

-0.0023 
+0.(3147 
-0.0057 
+0.0197 
-0.0080 
+0.0662 
-0.0369 


Mw selection 


15-25 
25-35 
35-65 
65-105 


+0.0030 

-0.0000 
+0.0018 
-0.0000 
+0.0068 
-0.0046 
+0.0000 
-0.0561 


+0.0000 

-0.0011 
+0.0010 
-0.0000 
+0.0034 
-0.0002 
+0.0032 
-0.0000 


+0.0000 

-0.0016 
+0.0034 
-0.0000 
+0.0000 
-0.0120 
+0.0560 
-0.0000 


+0.0001 

-0.0000 
+0.(J007 
-0.0000 
+0.0038 
-0.0000 
+0.0000 
-0.0598 


+0.0019 

-0.0000 
+0.0024 
-0.0000 
+0.0012 
-0.0040 
+0.0000 
-0.0149 


+0.0000 

-0.0007 
+0.(3(3(37 
-0.0000 
+0.0035 
-0.0000 
+0.0(313 
-0.(3210 


r background 


15-25 
25-35 
35-65 
65-105 


+0.0009 

-0.0009 
+0.0008 
-0.0010 
+0.0008 
-0.0012 
+0.0024 
-0.0025 


+0.00O1 

-0.0000 
+0.0000 
-0.0001 
+0.0001 
-0.0002 
+0.0003 
-0.0005 


+0.0014 

-0.0013 
+0.0023 
-0.0018 
+0.0018 
-0.0018 
+0.0000 
-0.0000 


+O.0001 

-0.0001 
+0.0002 
-0.0002 
+ 0.0002 
-0.0002 
+ 0.0000 
-0.0000 


+0.0011 

-0.0011 
+0.0013 
-0.0013 
+ 0.0011 
-0.0014 
+0.0015 
-0.0016 


+0.0001 

-0.0001 
+0.0001 
-0.0002 
+ 0.0001 
-0.0002 
+0.0002 
-0.0003 


Z background 


15-25 
25-35 
35-65 
65-105 


+0.0000 

-0.0009 
+0.0000 
-0.0010 
+0.0000 
-0.0012 
+0.0000 
-0.0025 


+0.0001 

-0.0000 
+0.0001 
-0.0000 
+0.0001 
-0.0000 
+0.0003 
-0.0000 


+0.0000 

-0.0017 
+0.0000 
-0.0037 
+0.0000 
-0.0029 
+0.0001 
-0.0000 


+0.0001 

-0.0000 
+0.0002 
-0.0000 
+ 0.0002 
-0.0000 
+0.0006 
-0.0007 


+0.0000 

-0.0012 
+0.0000 
-0.0019 
+ 0.0000 
-0.0017 
+0.0000 
-0.0016 


+0.0001 

-0.0000 
+0.00(31 
-0.0000 
+0.0(301 
-0.0000 
+0.0004 
-0.(3001 


QCD background 


15-25 
25-35 
35-65 
65-105 


+0.0140 

-0.0127 
+0.0092 
-0.0326 
+0.0018 
-0.0149 
+0.0134 
-0.0050 


+0.0009 

-0.0009 
+0.0013 
-0.0003 
+0.0010 
-0.0001 
+0.0002 
-0.0014 


+0.0329 

-0.0125 
+0.0242 
-0.0076 
+0.0062 
-0.0087 
+0.0243 
-0.0000 


+0.0010 

-0.0030 
+0.0009 
-0.0039 
+0.0013 
-0.0015 
+0.0039 
-0.0007 


+0.0206 

-0.0126 
+0.(3143 
-0.0242 
+0.0(300 
-0.0086 
+0.0174 
-0.0031 


+0.0009 

-0.0015 
+0.(3(312 
-0.0014 
+0.0004 
-0.0000 
+0.0006 
-0.(3001 


FMC signal event yield 


15-25 
25-35 
35-65 
65-105 


+0.0031 

-0.0065 
+0.0190 
-0.0000 
+0.0173 
-0.0000 
+0.0000 
-0.0540 


+0.0000 

-0.0023 
+0.0003 
-0.0016 
+0.0023 
-0.0033 
+0.0009 
-0.0052 


+0.0000 

-0.0004 
+0.0355 
-0.0000 
+0.0030 
-0.0022 
+0.0830 
-0.0000 


+0.0012 

-0.0032 
+0.0020 
-0.0017 
+0.0039 
-0.0000 
+ 0.0000 
-0.0557 


+0.0019 

-0.0043 
+0.0246 
-0.0000 
+ 0.0115 
-0.0000 
+ 0.0110 
-0.0267 


+0.0000 

-0.0016 
+0.00(38 
-0.0016 
+0.0028 
-0.0022 
+0.0000 
-0.0246 


W+7 


15-25 
25-35 
35-65 
65-105 


+0.0000 

-0.0060 
+0.0000 
-0.0064 
+0.0103 
-0.0000 
+0.0000 
-0.0704 


+0.0000 

-0.0031 
+0.0019 
-0.0000 
+0.0023 
-0.0000 
+0.0024 
-0.0000 


+0.0022 

-0.0000 
+0.0050 
-0.0000 
+0.0026 
-0.0000 
+0.0000 
-0.0060 


+0.0021 

-0.0000 
+0.0000 
-0.0039 
+0.0015 
-0.0000 
+0.0000 
-0.0166 


+0.0000 

-0.0031 
+0.0000 
-0.0025 
+0.0080 
-0.0000 
+ 0.0000 
-0.0467 


+0.0000 

-0.0015 
+0.0001 
-0.0000 
+0.0020 
-0.0000 
+0.000(3 
-0.(3049 


Aq and A4 variation 


15-25 
25-35 
35-65 
65-105 


+0.0296 

-0.2024 
+0.0691 
-0.2326 
+0.1053 
-0.1610 
+0.1189 
-0.0708 


+0.0074 

-0.0015 
+0.0155 
-0.0051 
+0.0218 
-0.0158 
+0.0150 
-0.0277 


+0.0323 

-0.1965 
+ 0.0553 
-0.1745 
+0.0917 
-0.1390 
+0.0963 
-0.0553 


+0.0080 

-0.0016 
+0.0240 
-0.0081 
+0.0170 
-0.0122 
+0.0167 
-0.0312 


+0.0305 

-0.2003 
+0.0644 
-0.2129 
+0.1013 
-0.1545 
+0.1106 
-0.0651 


+0.0076 

-0.0015 
+0.0181 
-0.0060 
+0.0202 
-0.0146 
+0.0157 
-0.(3290 


PDF variation 


15-25 
25-35 
35-65 
65-105 


+0.0000 

-0.0159 
+0.0152 
-0.0000 
+0.0000 
-0.0068 
+0.0022 
-0.0000 


+0.0013 

-0.0000 
+0.0000 
-0.0034 
+0.0011 
-0.0000 
+0.0134 
-0.0000 


+0.0000 

-0.0419 
+0.0000 
-0.0018 
+0.0759 
-0.0000 
+0.4001 
-0.0000 


+0.0017 

-0.0000 
+0.0000 
-0.0003 
+0.0000 
-0.0003 
+0.0000 
-0.0931 


+0.0000 

-0.0251 
+0.0094 
-0.0000 
+0.0178 
-0.0000 
+0.1486 
-0.(3000 


+0.0014 

-0.0000 
+ 0.0000 
-0.0024 
+0.0(306 
-0.0000 
+0.0000 
-0.0275 


Q 2 variation 


15-25 
25-35 
35-65 
65-105 


+0.1043 

-0.0000 
+0.1319 
-0.0000 
+0.1347 
-0.0000 
+0.0000 
-0.0959 


+0.0028 

-0.0000 
+0.0000 
-0.0104 
+0.0000 
-0.0267 
+0.0101 
-0.0000 


+0.1200 

-0.0000 
+0.0668 
-0.0000 
+0.0386 
-0.0000 
+0.0000 
-0.2251 


+0.0000 

-0.0121 
+ 0.0000 
-0.0076 
+0.0000 
-0.0180 
+0.0162 
-0.0000 


+0.1098 

-0.0000 
+0.1099 
-0.0000 
+0.1062 
-0.0000 
+0.0000 
-0.1434 


+0.0000 

-0.0018 
+0.0000 
-0.0096 
+0.000(3 
-0.0239 
+0.0124 
-0.(3000 
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TABLE XV: Total systematic uncertainties in the measure- 
ment of A 2 and A3, for electron and muon W+jet events and 
the combination of the electron and muon results. The sys- 
tematic uncertainties of Table XIV are combined in quadra- 
ture. 



Total systematic uncertainties 



p¥ (GeV) 


Electrons 


Muons 


Combination 


&A 2 


5A 3 


SA 2 


6A3 


SA 2 


6A3 


15-25 
25-35 
35-65 
65-105 


+U.1276 

-0.2209 
+0.1779 
-0.2622 
+0.1779 
-0.1699 
+0.2743 
-0.4229 


4-U.UUy2 

-0.0087 
+0.0191 
-0.0131 
+0.0305 
-0.0321 
+0.0883 
-0.0587 


+u.iayu 
-0.2220 

+0.1492 
-0.2098 
+0.1276 
-0.1643 
+ 0.5503 
-0.3941 


+U.U128 

-0.0140 
+0.0362 
-0.0208 
+0.0263 
-0.0243 
+0.0453 
-0.1629 


+U.13UY 

-0.2189 
+0.1641 
-0.2428 
+0.1530 
-0.1665 
+0.3297 
-0.3977 


+0.0094 

-0.0064 
+0.0239 
-0.0137 
+0.0289 
-0.0293 
+0.0692 
-0.0764 



TABLE XVI: The electron and muon CDF data event yields 
for inclusive VT+jet production, with statistical and system- 
atic uncertainties. The systematic uncertainties are due to 
By* and ?7j^ b cuts. 



Data event yields for inclusive W+jet production 


Pt (GeV) 


N e 




15-25 
25-35 
35-65 
65-105 


5166 ± 72 ±53g 
3601 ± 60 till 
3285 ± 57 tf 6 
624 ± 25 t\ 


2821 ± 53 tloi 
1869 ± 43 tH 
1880 ± 43 t\l 
371 ± 19 t° 



merit is reduced by more than 50%. Nevertheless, we 
prefer not to constrain the dataset in such a manner. 
Note that the event yield measurements do not affect the 
angular coefficient measurements, since we are interested 
only in the shapes of the distributions, in the latter case. 

X. FINAL RESULTS 

Combining the statistical and systematic uncertainties 
associated with the A2 and A3 measurement, we obtain 
our final results, presented in Tables XVII and XVIII. 
Figure 31 shows the measurement of A 2 and A 3 for the 
electron IF+jet data and Figure 32 shows the measure- 
ment of A2 and A3 for the muon PF+jet data. The com- 
bination of the electron and muon measurements of the 
two angular coefficients is presented in Figure 33. The 
Standard Model predictions for these angular coefficients, 
up to order a 2 s , are also presented. 

XI. CONCLUSIONS 

We have made the first measurement of the A2 and 
A3 angular coefficients of W boson production and de- 
cay, using the CDF Run la and Run lb electron and 
muon VF+jet data. Our datasets include at least one jet, 
satisfying the energy and pseudorapidity requirements. 
Due to finite statistical analyzing power of our VF+jet 
datasets and the characteristics of the W decay, only the 
measurement of A2 and A3 angular coefficients is sta- 
tistically significant, with the analysis of the azimuthal 



TABLE XVII: The measurement of the A 2 coefficient along 
with the statistical and systematic uncertainties for electron 
and muon W+jet events and the combination of the electron 
and muon results. The SM values up to order a 2 a are also 
included. 



Measurement of A 2 angular coefficient 




Pt (GeV) 


A 2 


Electrons 


15-25 
25-35 
35-65 
65-105 


0.02 ± 0.14 t u 2 i, 
0.14 ± 0.15 tlfe 
0.45 ± 0.13 t o\r 
1.24 ± 0.29 to.ll 


Muons 


15-25 
25-35 
35-65 
65-105 


0.14 ± 0.19+0^ 
0.55 ± 0.21 ti ll 
0.55 ± 0.20 to.ll 
0.98 ± 0.38lo;|9 


Combination 


15-25 
25-35 
35-65 
65-105 


0.06 ± 0.11^22 

0.28 ± 0.121°;^ 
0.48 ± 0.11 to.lr 
1.15 ± 0.23lo;io 


Standard Model 


15-25 
25-35 
35-65 
65-105 


0.09 
0.19 
0.35 
0.60 



TABLE XVIII: The measurement of the A3 coefficient along 
with the statistical and systematic uncertainties for electron 
and muon W+jet events and the combination of the electron 
and muon results. The SM values up to order o? s are also 
included. 



Measurement of A3 angular coefficient 




Pt (GeV) 


A 3 


Electrons 


15-25 
25-35 
35-65 
65-105 


0.03 ± 0.06 toli 
0.07 ± 0.06 toll 
0.13 ± 0.07 +o 03 
0.21 ± 0.15 till 


Muons 


15-25 
25-35 
35-65 
65-105 


0.03 ± 0.09+;!;'oi 
0.09 ± 0.091^2 
0.13 ± O.Wtlll 
0.33 ± 0.19 ±0+6 


Combination 


15-25 
25-35 
35-65 
65-105 


0.03 ± 0.05 +0.01 
0.08 ± 0.05l^i 
0.13 ± 0.06 ttol 
0.26 ± 0.12 till 


Standard Model 


15-25 
25-35 
35-65 
65-105 


0.02 
0.04 
0.08 
0.16 



angle of the charged lepton in the W rest-frame. The 
A and A 4 coefficients are preferably measurable with a 
polar angle analysis, while A\ and the next-to-leading 
order coefficients - A5, Aq, and A7 - are not measurable, 
with any meaningful statistical significance, with Run I 
VF+jet data. 

At leading order, the A 2 and A 3 angular coefficients 
fully describe the azimuthal W angular distribution in 
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FIG. 31: Measured A2 and A3 using electron W+jet events. 
The total (outer) and statistical (inner) uncertainties are 
shown along with the Standard Model 1-loop prediction up 
to order a 2 (dashed line). 



FIG. 32: Measured A2 and A3 using muon VF+jet events. The 
total (outer) and statistical (inner) uncertainties are shown 
along with the Standard Model 1-loop prediction up to order 
a 2 (dashed line). 



the Collins-Soper W rest-frame. These angular coeffi- 
cients are also part of the total W differential cross sec- 
tion, and can be expressed as ratios of the correspond- 
ing helicity cross sections of the W to its total unpolar- 
ized cross section. This measurement tests the Standard 
Model prediction for W polarization, and the associated 
QCD corrections present in the production of W bosons 
at high transverse momenta. We observe good agree- 
ment with the Standard Model prediction up to order o? s 
in QCD. 
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